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Abstract 

Characteristic features of transverse transport along the a* axis in the NbSe3 charge-density- wave 
conductor are studied. At low temperatures, the I-V characteristics of both layered structures and 
NbSe3-NbSe3 point contacts exhibit a strong peak of dynamic conductivity at zero bias voltage. 
In addition, the I-V characteristics of layered structures exhibit a series of peaks that occur at 
voltages equal to multiples of the double Peierls gap. The conductivity behavior observed in the 
experiment resembles that reported for the interlayer tunneling in Bi-2212 high-T c superconductors. 
The conductivity peak at zero bias is explained using the model of almost coherent interlayer 
tunneling of the charge carriers that are not condensed in the charge density wave. 

PACS numbers: 71.45.Lr; 72.15.Nj; 74.50. +r 
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It is well known that the crystal structure of BSCCO high-T c superconductors consists of 
atomically thin superconducting cuprate layers spatially separated by atomically thin insu- 
lating layers of BiO and SrO. The corresponding spatial modulation of the superconducting 
order parameter in the direction across the layers (along the c axis) leads to a discrete de- 
scription of the transverse transport (the Lawrence-Doniach modelB) with the neighboring 
superconducting layers being coupled through Josephson tunneling junctions. The validity 
of this approach is confirmed by the results of the experiments on natural Bi-2212 layered 
structures of small lateral size!!. Currently, the study of the interlayer tunneling of Cooper 
pairs and quasiparticlesB! is one of the new original methods for investigating high-T c su- 
perconductors. 

In this paper, we study the interlayer tunneling in a layered system with an electron 
condensate of a different type, namely, in a charge-density-wave (CDW) conductor. The 
material chosen for our experiments is NbSe3. This compound is characterized by two 
Peierls transitions, which occur at the temperatures T p \ =145 K and T P 2 =59 K. In the 
low-temperature Peierls state, the Fermi surface retains some regions where the nesting 
conditions are not satisfied (the pockets) and, hence, the Peierls gap is absent. Therefore, 
NbSe3 does not undergo transition to the insulating state and retains metallic properties 
down to the lowest temperatures^. From the analysis of both the crystal structure of NbSe3 
and the anisotropy of its properties, it follows that this material can be classed with quasi- 
two-dimensional layered compounds. In fact, its conductivity anisotropy in the (b — c) plane 
is determined by the chain conductivity along the b axis and is estimated as <Jb/<J c ~ 10, 
whereas the conductivity ratio o"b/cr a » is determined by the layered character of the structure 
and reaches the values ~ 10 4 at low temperatures^'!. Figure la shows the crystal structure 
of NbSe3 in the (a — c) plane. One can see that, in this material, the layered structure 
is formed as a result of a pairwise coordination of selenium prisms with the predominant 
orientation of their bases in the (b — c) plane. In Fig. lb, the shaded areas indicate the 
elementary conducting layers in which the prisms are rotated and shifted with their edges 
toward each other. In these layers, the distances between the niobium chains are relatively 
small, whereas the neighboring conducting layers are separated by an insulating layer formed 
as a double barrier by the bases of the selenium prisms. 

This type of layered structure in combination with the conductivity anisotropy offers the 
possibility for the CDW condensation in the elementary conducting layers spatially separated 
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by atomically thin insulating layers. In this case, as in layered high-T c superconductors, one 
can expect that the CDW order parameter will be modulated along the a* axis, and the 
transport across the layers will be determined by the intrinsic interlayer tunneling between 
the elementary layers with the CDW. 

To verify these speculations, we carried out an experimental study of the transport across 
the layers in NbSe3 in the condensed CDW state. 

The samples used in the experiment were layered structures with a small area through 
which the current flows across the layers, S = 2 x 2 /im 2 , and with the number of elementary 
layers ~30 (the overlap junctions). They were made from a thin NbSe3 single crystal by 
the focused ion beam processing technique developed for fabricating similar structures from 
Bi-2212!. Schematically, the structure under investigation is shown in the inset of Fig. 2. 
In addition, we studied the characteristics of NbSe3-NbSe3 point contacts oriented in the 
direction of the a* axis. The configuration of such a contact is shown in the inset of Fig. 
3. The contact was formed directly at low temperature by bringing two NbSe3 whiskers 
together with high accuracy. In both cases, we performed four-terminal measurements of 
the I-V characteristics and their derivatives at low temperatures, i.e., below the second 
Peierls transition T p2 = 59 K. 

The main results of the experiment are shown in Fig. 2, which, for one of the virtually 
perfect structures (No. 3), displays the dynamic conductivity along the a* axis versus the 
bias voltage V at different temperatures from 59.5 to 4.2 K. It can be seen that, for T < T P 2, 
the differential I-V characteristics are of a tunneling character. When the temperature 
decreases below ~35 K, the I-V characteristics begin to exhibit a conductivity peak at zero 
bias voltage. As the temperature decreases further, this peak becomes a dominant feature. 
The peak amplitude reaches saturation when the temperature decreases below 6-8 K, and at 
T = 4.2 K, it is almost 20 times as great as the conductivity observed at high bias voltages. 
Note that this anomaly cannot be attributed to Joule heating. The estimate of the sample 
heating for a typical value of heat transfer to helium yields a conductivity decrease of less 
than 10% for the whole range of measured voltages. In addition, the I-V characteristics 
show a clearly defined set of conductivity peaks that are symmetric about zero voltage and 
at low temperature correspond to \V\ =50, 100, and 150 mV, i.e., to |V| = nV , where V Q = 
50 mV and n= 1,2, 3. As the temperature increases above 25 K, these peaks move to lower 
energies, and, at the temperature T ^59 K corresponding to the second Peierls transition 
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for NbSe3, the values of V n vanish. 

The picture described above was observed only for perfect structures. The presence 
of structure defects, such as a twin boundary, leads to a considerable reduction of the 
conductivity peak at zero bias and to the appearance of a peak at V — 25-27 mV ~ Vq/2 
(curve 2 in Fig. 3). Qualitatively similar dependences were observed for some of the NbSe3- 
NbSe3 point contacts. Their I-V characteristics also exhibit a conductivity peak at zero 
bias with an amplitude approximately equal to that observed for defect overlap junctions, 
as well as a conductivity peak at V ~ 25 mV (curve 1 in Fig. 3). Figure 4 presents the 
temperature dependences of the normalized dynamic conductivity at zero bias voltage for 
all types of samples studied in the experiment. One can see that, for a perfect overlap 
junction, the amplitude of the zero bias anomaly is more than three times greater than 
the amplitudes observed for the defect layered structures and point contacts. Note that 
the voltage value Vq = 50 mV is close to twice the value of the low-temperature energy 
gap, 2A p2 /e, for NbSe30'0, while the dependence V n (T)/V n (0) obtained for both overlap 
junctions and point contacts (Fig. 5) agrees well with the temperature dependence of the 
energy gap predicted by the BCS theory (the dashed line in Fig. 5). This result suggests 
that the conductivity features observed in the experiment are governed by a gap mechanism. 

It is important to note that all effects described above were observed exclusively for the 
transport across the layers (along the a* axis). The extra experiments performed by us on 
specially fabricated bridges and point contacts oriented along the c axis (with the transport 
across the chains in the layer plane) showed no conductivity peak at V = 0. 

Let us first analyze the results obtained for perfect overlap junctions. The most prominent 
feature of the I-V characteristics of these structures is the strong conductivity peak at 
zero bias voltage. In addition, the characteristics exhibit a periodic sequence of peaks 
at \V\ = uVq, which resembles the series of quasiparticle branches observed in the I-V 
characteristics of Bi-2212 layered structures when the measurements are performed across 
the layers!. As was noted above, the value Vq = 50 mV is close to twice the value of the 
low temperature CDW gap in NbSe3, and the temperature dependence Vo(T) follows the 
prediction of the BCS theory. However, in contrast to the case of a superconductor, the 
conductivity observed for NbSe3 at zero bias voltage is finite. This fact suggests that the 
interlayer conduction at zero bias does not originate from the collective CDW contribution 
of the Josephson tunneling type. 
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The conductivity peak observed at zero bias cannot be explained by the regular (inco- 
herent) single-particle tunneling as well. If this were the case, the conductivity peak could 
be attributed to the energy dependence of the density of states of electrons that are not 
condensed in the CDW. Then, the conductivity feature under discussion should also be ob- 
served for N-I-CDW tunnel junctions. However, no such effects were revealed by the detailed 
experimental studies of this kind of tunnel junctions oriented along the a* axi&BS. 

We believe that the conductivity peak at zero bias is caused by the coherent interlayer 
tunneling of the charge carriers localized in the pockets of the Fermi surface, and below we 
thoroughly substantiate this statement. 

In the general case, the tunnel current between two nonsuperconducting layers, a and b, 
is described by the expression^!: 

1 = ^JT E MP> / + duJ ( tanh Jf ~ tanh ^T~) ImG - (p ' ^ G b(^ w - eV), 

(1) 

where G R is the retarded Green function of an electron with the momentum p in the 
layer a, lmG R (p,u) is the spectral density of the Green function, t ab (p,q) is the matrix 
element characterizing the tunneling from the state with the momentum p in the layer a to 
the state with the momentum q in the layer b, V is the voltage between the adjacent layers, 
and T is the temperature. 

Note that the contribution of the collective (moving) CDW modei to the interlayer 
tunneling is absent, because the tunnel current is determined by the Green functions in 
different layers, whereas the collective mode propagates within the layers. 

Below, we consider the coherent tunneling, when the electron momentum in the plane 
does not change under the tunneling: p = q. Then, we have 

/= dp\t ab (p)\ 2 J + du (tanh^-tanh^^) hnG R (p, oj)G r (p, oj - e V) (2) 

Taking into account that the scattering in a layer is determined by the collision frequency 
v, we obtain 



ImG^p, oj) = 7 2^ ( 3 ) 

7r(u; — e{py + j 2 ) 
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where e(p) is the electron spectrum and 7 = fiv. We can use this expression with 7 = 
lsc + line, which takes into account the change in the momentum because of the scattering 
within the layers and the change due to the tunneling ^ inc = (e(p) — e(q)). In both cases, 
we take into account the energy uncertainty for the state characterized by the momentum p. 
We define the tunneling as almost coherent when ^ inc and 7 are small relative to other energy 
parameters of the electron system. In our case, such parameters are the Peierls gap and the 
width of the electron band in the pockets. Replacing £ a &(p) by the momentum- independent 
quantity t and integrating with respect to p, we obtain 

^ / rfp|U(p)| 2 ImGf(p,^)ImGf(p,^- eV) = fffi+^I ( 4 ) 

where N(0) is the density of states of electrons in the pockets. Finally, for the interlayer 
current when eV < 2 A, we obtain the expression 

ey 7 iV(0)|t| 2 

1 ' 27r 3 (e 2 ^ 2 + 4 7 2 ) 1 ' 

and the dynamic conductivity has the form 

^g)=47 2 4 ^ (6) 
<r(0) 1 (e 2 U 2 + 4 7 2 ) 2 1 ' 

Note that the temperature dependence is present in 7(T) only. One can see that the 
dynamic conductivity has a peak of width m 7 at V = and becomes negative and unstable 
when eV > 27. For a thirty-layer structure, the experimentally observed peak width ~10 mV 
corresponds to 7 ~ 0.3 meV. For comparison, we estimate the parameter ^ sc characterizing 
the scattering within the layers. Using the known mobility data fi = e/v sc m* ~ 4 x 10 4 
cm where m* = 0.24m e li3, we obtain 7 SC w 0.13 meV; i.e., the changes in the electron 

momentum because of tunneling are approximately equal to changes caused by the scattering 
within the layers. 

Now, we explain why the interlayer current decreases with increasing voltage V. The 
electron tunneling between the layers must obey the energy conservation law; i.e., e(p) = 
e(q) — eV to within 27. For coherent tunneling, we have e(p) = e(q), and for eV 7, 
tunneling is possible and we have the conventional Ohm's law. When eV > 27, tunneling 
is impossible up to the voltage V reaching 2A/e. In this case, electrons condensed in the 
CDW begin to contribute to the interlayer current in the form of a regular tunneling of 
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CDW quasiparticles through the double Peierls gap 2 A. Thus, the interlayer current can be 
realized by means of only one of the two aforementioned mechanisms. In an actual multilayer 
structure, because of the geometric nonuniformity of individual layers (the areas of the layers 
are somewhat different), the voltage value V ~ 27/e will not be reached simultaneously for 
different individual tunneling junctions. When V > 27/e, some of the junctions can be in 
the coherent tunneling regime whereas other junctions can be in the regime of single-particle 
tunneling through the gap. As a result, a sequence of conductivity peaks must appear in 
the I-V characteristic of the compound under study at the voltages V = 2nA p /e, where n 
= 1,2,.... 

In the presence of defects in an overlap junction or in the case of a point contact, the 
incoherence of tunneling is enhanced, although in the best structures the tunneling remains 
almost coherent. Therefore, the I-V characteristics of these structures, along with the 
conductivity peak at zero bias, contain an additional feature at V = A p /e, which is related 
to the single-particle tunneling of the N-I-CDW type. 

Thus, the results of this study show that the interlayer tunneling in natural layered 
structures obtained from NbSe3 represents an independent efficient method for investigating 
the Peierls state in this compound. The tunneling conductivity peak observed at zero 
voltage can be self-consistently explained by the almost coherent interlayer tunneling of 
charge carriers that are not condensed in the CDW and are localized in the Fermi surface 
pockets not covered by the gap. The series of equidistant peaks in the I-V characteristic 
can be explained by the quasiparticle tunneling through the CDW gap under sequential 
transitions of individual tunneling junctions to the resistive state. Note that the conductivity 
peak observed at zero bias is a unique manifestation of coherent single-particle transport in 
solids. 

We are grateful to S.A. Brazovskii, V.A. Volkov, and V.M. Yakovenko for useful dis- 
cussions. This work was supported by the Russian Foundation for Basic Research (project 
nos. 99-02-17364, 01-02-16321, and 00-02-22000 CNRS), the state program Physics of Solid 
Nanostructures (project no. 97-1052), the NWO Russian-Netherlands Project, and the Los 



7 



Alamos National Laboratory with the support of US DOE. 



1 W. E. Lawrence and S. Doniach, in Proceedings of the 12th International Conference on Low 
Temperature Physics, Kyoto, Japan, 1971, Ed. by E. Kanda, p. 361. 

2 R. Kleiner and P. Miiller, Phys. Rev. B 49, 1327 (1994). 

3 R. Kleiner, P. Miiller, H. Kohlstedt, et al., Phys. Rev. B 50, 3942 (1994). 

4 K. Lanabe, Y. Hidaka, S. Karimoto, and M. Suzuki, Phys. Rev. B 53, 9348 (1996). 

5 Yu. I. Latyshev, L. Yamashita, L. N. Bulaevskii, et al., Phys. Rev. Lett. 82, 5345 (1999). 

6 G. Griiner, Density Waves in Solids (Addison - Wesley, Reading, Massachusetts, 1994) 

7 N. P. Ong and J. W. Brill, Phys. Rev. B 18, 5265 (1978). 

8 Yu. I. Latyshev, P. Monceau, O. Laborde, et al., J. Phys. IV 9, 165 (1999). 

9 S.-J. Kim, Yu. I. Latyshev, and T. Yamashita, Appl. Phys. Lett. 74, 1156 (1999). 

T. Ekino and J. Akimitsu, Physica B 194-196, 1221 (1994). 

1 A. Fournel, J. P. Sorbier, M. Konczykovski, and P. Monceau, Phys. Rev. Lett. 57, 2199 (1986). 

2 I. O. Kulik and I. K. Yanson, Josephson Effect in Superconducting Tunnel Structures (Nauka, 
Moscow, 1970). 

3 P. A. Lee, T. M. Rice, and P. W. Anderson, Solid State Commun. 14, 703 (1974). 

4 N. P. Ong, Phys. Rev. B 18, 5272 (1978). 

5 R. V. Coleman, M. P. Everson, Hao-An Lu, et al., Phys. Rev. B 41, 460 (1990). 



8 



a) 



y=0 b/2 



Se 
Nb 



6) 



FIG. 1: Schematic representation of the NbSe3 structure in the (a — c) plane; (b) the same 
structure with the conducting planes indicated by shading 
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FIG. 2: Dependences of dl/dV on the voltage V for an overlap junction (sample No. 3) at different 
temperatures: 59.5, 55.6, 52.1, 48.0, 43.7, 40.1, 37.0, 34.3, 31.0, 28.0, 25.1, 22.6, 19.1, 15.7, 12.8, 
8.0, and 4.2 K. The dynamic conductivity scale corresponds to the curve at T = 59.5 K, and other 
curves are shifted upwards for clarity. The inset shows the sample configuration. 
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FIG. 3: Dependences dI/dV(V) (1) for a NbSe3-NbSe3 point contact and (2) for an overlap junction 
(sample No. 2) at T = 4.2 K. The inset shows the point contact configuration. 




FIG. 4: Temperature dependences of the dynamic conductivity at zero bias voltage for two overlap 
junctions: sample No. 3 (full circles) and sample No. 1 (full squares), and for a NbSe,3-NbSe3 
point contact (empty circles). The conductivity values are normalized to the value of dl/dV (T = 
62 K). 
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FIG. 5: Temperature dependence of the second Peierls energy gap for an overlap junction (sample 
No. 3, full circles) and a NbSe3-NbSe3 point contact (empty circles). 
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